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-* Introduction 
¢ Schematic of tube-on-tank liquefaction concept 


¢ Modeling of tube-on-tank configuration 
LY Modeling case: 


v Incoming gas with temperature of 273 K (warm case, 
baseline) and 100 K (cold case) 


Y Tank fill level: 0% and 95% 
LJ Modeling approach: 
v 2D axisymmetric CFD model in ANSYS Fluent 


> Investigate the mixing of incoming gaseous O2 with the 
fluid inside the tank 


¥Y 1D thermal model in Matlab 
> Understand how to set the BCs in the thermal model 
¥Y 3D thermal model in MSC Patran/Pthermal 


> Investigate the thermal gradient near the top of the tank 
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Introduction 


¢ The in-situ production of propellants for Mars 
missions will utilize Mars atmospheric carbon dioxide 
(CO.) to produce oxygen. 
¢ The oxygen is then cooled, liquefied, and stored to be 
available for Mars ascent propulsion system, which 
could be up to 2 years after liquefaction starts. 
¢ Recent investigations have demonstrated the 
feasibility of using high-efficiency reverse turbo- 
Brayton-cycle cryocoolers to: 1*'stage tank 
- Cool the oxygen gas 
¢ Liquefy the oxygen gas ion 


2d stage 
i radiator 


¢ Achieve zero boil-off 
¢ Control the pressure of we 
oxygen within a tank Helium tank 


& dic ||| ; 


2"d stage 
tank 
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NASA 


in-situ Production — Liquefaction - Storage 


DRY GOX 
at 273 K and 1 atm 


al 


p 


MARS ENVIRONMENT 


Reference: 

1. “Mars Ascent Vehicle Design for Human Exploration”, Tara Polsgrove, AIAA SPACE 
2015. 

2. “MAV Deep Dive: ISRU to MAV Propulsion Interface, Update on LOX Production, 
Liquefaction and Transfer v2.0”, Bill Studak, Aug. 15 2015. 
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Cryocooler 
system _....... 


Turbine with 
Alternator 


Cooling tubes 


Compressor/Motor 
Assembly 


: eo ry ae 
oy Sage ee eae = _ 
t 
wren tegen — a 3 rots bot ety 
eet =~ Ske ae tea oo 


Radiator with 
Heat Pipes 


LS She oe ae ae 
ae 
fe mcenr 


. Oo See 
\ ~~ 


Tank Support Struts 


The gaseous neon circulating in the cryocooler system is maintained 
slightly below liquid oxygen saturation tem perature and is routed 
through a network of cooling tubes epoxied to the tank wall. 

The oxygen gas produced from the in-situ production process is 
introduced into the chilled tank. 
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CFD Model in ANSYS Fluent 


¢ 2D axisymmetric, transient analysis 
LJ Multiphase model: Mixture/slip 
velocity/implicit body force 


ae 


LI Turbulence model: shear stress transport ee 
(SST) k-w (2 eqns) ; siiteriisee? 

- No conjugate heat transfer (Tank wall g sunioat 
and neon tubes are not modeled) sutcereseas’ 

U Simplify Fluent CFD model to save suena 


a 
ra; 


computational time 
UI Define tank wall boundary condition 
(constant T at 90 K or heat flux at - 12 W/m? 
= - 243.6 W/20.3 m2 based on lift of 
cryocooler) cr 
UU Investigate uncertainty of decoupling neon 
cooling tube and tank wall 


Guu een. 
Ky 
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¢ Fluent model results will be shown for 
¢ Filllevel: 0% and 95% 
¢ Incoming warm GOX at the mass flow rate of 
2.2 kg/hr 
¢ Incoming pre-chilling GOX at the mass flow 
rate of 2.2 kg/hr 
¢ Wall boundary conditions: 
(a) constant tank wall temperature 
(b) constant tank wall heat flux 
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Temperature contour of mixture of GOX and LOX 
o Incoming gas: 273 K 
(a) wall temperature fixed at 90 K 
(b) wall heat flux fixed at -12 W/m? 


Le 


Contours of Static Temperature (mixture) (k) (Time=2.3095e+03) Jun 14, 2016 Contours of Static Temperature (mixture) (k) (Time=2.4000e+03) Jun 15, 2016 
ANSYS Fluent Release 16.2 (axi, dp, pbns, mixture, sstkw, transient) ANSYS Fluent Release 16.2 (axi, dp, pbns, mixture, sstkw, transient) 


¢ Thewarm gaseous O, chills down within smaller volume with a cold 
wall as shown in case (a). 
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Time history of the mass of LOX: 
o Incoming gas: 273 K and 100 K 
(a) wall tem perature fixed at 90 K 
b) wall heat flux fixed at -1 2 W/m2 


Tgas = 273K 
12 —— Tgas = 100 K 


mass of LOX (kg) 
S 
o 


0 5 10 15 20 25 30 35 40 


0 5 10 15 20 25 30 35 40 
(min) 


time 


¢ TheLOX mass produced inside the tank at t = 40 minutesis 
¢ Forincoming gas of 273 K: 
¢ 1.48 kg in case (a), 0.55 kg in case (b), a ratio of 2.7. 
¢ Forincoming gas of 100 K: 
¢ 1.52 kg in case (a), 0.95 in case (b), a ratio of 1.6. 
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* centerline 
. 


( a) Wall pool, 


2.00e+02 


temperature wy 


fixed at 90 K * nl 


1.00e+02 — 


8.00e+01 


(a) 


ch , 


0.75 1 1.25 15 1.75 
Position (m) 


* centerline 


2.80e+02 

2.600+02 \ 
4 . 
2.40e+02 — z 


2.20e+02 4 


(b) Wall heat —— 


Temperature 
a iienene: 


flux fixed at -12 ‘ 


1400+02 


1.20e+02 + 


W/m2 = 


Static Temperature (mixture) (Time=2.3095e-: 


+03) Jun 16, 2016 
ANSYS Fluent Release 16.2 (axi, dp, pbns, mixture, sstkw, transient) 


II 


oO 0.25 05 0.75 1 1.25 15 1.75 2 2.25 25 


Position (m) 


8/10/2017 


Static Temperature (mixture) (Time=2.4000e+03) 


Jun 16, 2016 
ANSYS Fluent Release 16.2 (axi, dp, pbns, mixture, sstkw, transient) 
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emperature contour of the a t= 20 mins 


2.55e+02 
. 2.466+02 
mixture of GOX and LOX for | =: 
2.27e+02 
1. 1. 2.18e+02 
incoming gas at 2/3 K — 
Q Q 2.00e+02 
1.916+02 
1.81e+02 
1.72e+02 
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Tape t= 0 min, Initial T ae 
Bete 3 1 45e+02 
9.90e+01 . . 41.36e+02 
9.85e+01 inside the tank Teter Oe 
SaGSEOT 1.17e+02 
er 1.08e+02 
SSIS507 9.91e+01 
9.00e+01 0 1@n) 
9.65e+01 SSS 
9.60e+01 
9.55e+01 Contours of Static Temperature (mixture) (k) (Time=1.2000e+03) Feb 15, 2017 
9.506e+01 ANSYS Fluent Release 17.1 (axi, dp, pbns, mixture, sstkw, transient) 
9.45e+01 
9.40e+01 
9.35e+01 
9.30e+01 
9.25e+01 
9.20e+01 
9.15e+01 
9.10e+01 
9.05e+01 
9 00e+01 9 ————! §"? 
x 
Contours of Static Temperature (mixture) (k) (Time=0.0000e+00) Feb 16, 2017 ° 
ANSYS Fluent Release 17.1 (axi, dp, pbns, mixture, sstkw, transient) S 
0) 
w 
E 


Tank wall boundary condition 
doesn’t change the liquefaction 
rate for 95% fill level case. 


time (mins) 
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LN Observation from Fluent Results ASA 


¢ Fluent model results show the mixing of the warm 
incoming GOX with the gas inside the tank. 


¢ Fluent results provide temperature distribution of 
incoming warm gas. 

¢ Tank wall boundary conditions show significant 
difference of liquefaction rate for 0% filllevel, but very 
little difference for 95% fill level. 

¢ The entire picture of heat transfer from neon gas to the 
tank wall then fluids is not shown in Fluent analysis. It 
will be interesting to know temperature changes of the 
neon fluid along the tube and the temperature gradient 
near the top of the tank. 

¢ 1D thermal circuit is built to understand more of the 
tube-on-tank configuration. 
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ALYSIS WORKSHOP 


R,s the tank 


r, neon,2 rT. Neon,3 T neon,4 , neon,5 


Conduction resistance between the wall nodes along the axial/circumferential directions 
Convection resistance between the cooling tube wall and neon fluid 

Convection resistance between the tank wall and gaseous O2 

Contact resistance between the cooling tube and tank wall 

Twat and To», distribution are needed to specify as BC 

Inlet temperature of neon gas and mass flow rate need to be defined 
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”)% 1D Model Results (I) NASA 


° Twan= (Teast sy)/2.0 at the top is used, T,, is the saturated vapor temperature 

¢ Neon gas inlet temp is assumed to be 80 K 

- Estimate the tank surface area A needed to cool the warm gas (7,,,) to the 
saturated temperature UuSINg Mgo:0,(T gas-Tsv) = 1 A (Tgas-T wan), then compute 
the tank height (= 0.42 m) based on A, which is at 94% fill level assuming h = 
0.5 W/m2-k 


——®— nodes attached tube 
—@— nodes between tubes 
—@— Neon 

—@— Gas 02 


—®— nodes attached tube 
—@— nodes between tubes 
——@— Neon 

—@— Gas 02 


Tyas = 100 K 
Tyas = 273 K ) ) 


Tank Wall T (K) 


0 0.2 0.4 06 0.8 1 12 14 16 18 2 02 04 06 08 4 12 14 16 18 2 
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1D Model Results (Il) 


Twain = Tas, Same as the incoming 7,,.. 
- Inside the tank, assume the temperature of T,,.= 7, - 


—®— nodes attached tube 
—©— nodes between tubes 
—@— Neon 

—@— Gas 02 


—®— nodes attached tube 
—@— nodes between tubes 
—@— Neon 

—@— Gas 02 


Tgas = 100 K 


Tgas = 273 K 


Tank Wall T (K) 
Tank Wall T (K) 


01 O02 O03 04 O5 O06 O7 O8 O89 1 0 01 Oo2 O03 O04 O5 O06 O7 O08 09 1 
tank height (m) (from top to the bottom) tank height (m) (from top to the bottom) 
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| 4 N ummary Of 1D Tube-on-tank Model Results ASA 


o There are uncertainties on how to define the incoming 
GOX temperature distribution inside the tank and the 
tank wall temperature near the top of the tank. 

o 1D model can not accurately show the gradient since 
the mesh size is limited. 

o 1D thermal circuit model shows the major BCs and 
assumptions that need to be considered for the 
modeling. 

o 3D tube-on-tank model in MSC Patran/pthermal is 
built to investigate the temperature gradient on the 
top of the tank. 
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Neon Tube 


LOX 
(95% 
fill- 
level) 


8/10/2017 


Steady-state analysis 
Geometry: 60° wedge + 


of the MAV tank (6 
cooling tubes) 

FEM mesh for large 
temperature gradient 
on the top of the tank 
Conduction is 
modeled for both 
GOX and LOX 
Convection is not 
modeled, phase 
change is not 
modeled 
Temperature of 
incoming GOX from 
Fluent model is used 
as BC 


Raat 


Les 


FEM mesh 
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/ XN MSC Patran/Pthermal Tube-on-tank Model Results (I) ASA 


Apply the Fluent model results of the incoming GOX temperature 

along the center line of the tank 

¢ Specify the tank wall temperature at the top equal to the incoming 
hot gas temperature 

¢ Specify the neon gas inlet temperature and mass flow rate 


Patran 2012.2 64-Bit 16-May-17 15:07:17 2.73+002 
Fringe: TIME: 0.0000000000D+00 SECONDS S.S.ITERATIONS: 341, nr000.nrf05, Temperature. ( deg-K ). . (NON-LAYERED) 2 60+009 


2.48+002 
2.35+002 
2.23+002 
2.10+002 
1.98+002 
1.85+002 
1.73+002 
1.60+002 
1.48+002 
1.35+002 


1.23+002 


1.10+002 


9.75+001 
8.50+001 
default_Fringe : 


. rs Max 2.93+002 @Nd 14584 
Soles Min 8.50+001 @Nd 1 
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| Wall tem perature distribution along the Xloc (height of the 


tank) with and without neon cooling (the worst case) 


280 
© With Neon cooling 
250 + Without Neon Cooling 
240 
220 
2 
o 200 
=J 
5 Xloc = 6.876” is at the top of the tank 
E 160 


100 batde 


"9000006006 @ 


: ¢ " 12 14 16 18 20 
~98% filllevel = «(rc 
¢ Results show the temperature near 


the top of the tank cools to 90 K within 
a short distance 
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Xloc 
(inch) 


6.876 
6.911 
6.993 
7.124 
7.304 
7.535 
7.817 
8.151 
8.539 
8.980 
9.476 
10.028 
10.635 
11.299 
12.020 
12.802 
13.647 
14.556 
15.534 
16.580 
17.699 
18.892 


273.000 
134.346 
127.414 
114.595 
105.333 
105.297 
100.350 
97.158 
95.223 
93.760 
92.786 
92.056 
91.411 
90.967 
90.657 
90.431 
90.272 
90.172 
90.102 
90.049 
90.017 
90.000 


(with Neon) (No Neon) 


273.000 
157.024 
124.497 
110.559 
104.443 
99.997 
97.055 
95.246 
93.703 
92.574 
91.708 
91.214 
90.858 
90.598 
90.377 
90.233 
90.132 
90.070 
90.032 
90.011 
90.003 
90.000 


20 


7 4% Observations From 3D Patran Thermal Model Results (& aus 


¢ Patran/pthermal results show a clear picture of the 

temperature gradient near the top of the tank due to 
cold neon and incoming hot GOX. 

¢ The tank wall temperature drops to 90 K from 273 K 
within a short distance, that is above the 97-98% fill 
level, even for the case of no neon cooling. 

¢ Based on three model results, we can conclude that 
liquefying the warm GOX without pre-chillis feasible 
and no major concern near the top of the tank for the 
thermal gradient. 

¢ The liquefaction rate over long time period (42+ days) 
was investigated using a separate thermal model in 
Thermal Desktop/Sinda-Fluint. 
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